UNCLASSIFIED 

An  <06  371 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  sped* 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 

Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  In  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  Is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  In  any  way  be  related 
thereto. 


<f)  2#/  MO 

( /{^TECHNICAL  |k)TE  No.  19 


CpJ-J'P 


(7'0*& 

(2)09  Ul 

'  f>  . 

<B>v* 

(/G-H)a//4 

(io)  u_ 

(Jj)  ^  4 


/a 


M.  Gourdin,^ 


Facul-^e  des  Scjjmcea^.Qrsay.  France. 
Faculte  des  Sciences,  Bordeaux,  France, 


J.  Tran  Thanh  Van  . 

Faculty  des  Sciences,  Orsay,  France. 


The  research  reported  in  this  document  has  been  sponsored  in  part  by  the 
’Air  Research  and  Development  Command,  United  States  Air  Force',  under  contract 
AF  61  (052)  -  175,  through  the  European  Office,  ARDC. 


Hb _ 

I  DA 


CampMg- 


SUMMARY 

•me  research  on  the  deuteron  relativistic  wave  function- - 

.  v.  A 


liinm  iijiyl  hwi^Bethe  and  Salpeter'srelativistic  wave  equation^to  'the  bound  system 
of  two  particles  with  spin  g  .  Af^er^some  very  slight  approximations,  the 
solution  of  the  infinite  system  of  coupled  integral  equations  ifcaeTesias 


the  explicit  structure  of  the  wave  matrices  corresponding  to  the  S  and  D  states^ 
tc  (<~f  ,  The  percentage  of  D  state  obtained  is  in  good  agreement  with  the  experimental 

result. 
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I.  INTRODUCTION 


We  have  already  studied  the  problem  of  the  determination  of  a  relativis¬ 
tic  deutercm  wave  function  ^  .  In  this  paper,  we  investigate  the  case  of  two 

(2) 

Dirac  nucleons  with  the  Bethe  and  Salpeter  equation  formalism  .  The  expansion 

of  the  wave  function  in  hyperspherlcal  harmonics  has  been  first  introduced  in 

(3) 

the  nucleon-nucleon  scattering  problem  ,  It  was  used  in  our  first  paper  for 
a  spinless  deuteron.  We  employ  here  the  PS  (PS)  meson  theory  with  an  interaction 
invariant  under  rotations. 

We  obtain  an  infinite  coupled  system  of  radial  equations  for  the 
16  components  of  the  wave  function.  In  the  first  section,  we  develop  some 
calculations  in  order  to  obtain  the  equations  corresponding  to  the  J  =  1  bound 
neutron-proton  state.  Because  of  parity  conservation,  we  can  use  selection 
rules  to  simplify  the  equations* 

In  the  second  section,  we  analyze  the  structure  of  the  integral 
equations.  There  appear  three  quantum  numbers  L,  m,  n  due  to  the  expansion 
in  hyperspherlcal  harmonics,  of  the  three  angles  0,  iP  ,  P,  corresponding  to  the 
four-vector  p.  The  first  one  L  is  the  orbital  angular  momentum  ;  the  second  m, 

*V\ 

its  projection  over  Oz.  The  last  one  n  has  no  direct  physical  interpretation  ; 
its  parity  is  conserved  and  it  is  responsible  for  the  coupling  between  integral 
equations.  But  it  is  very  useful  to  note  that  this  coupling  is  extremely  weak  } 
practically  the  lowest  value  for  n  is  important  and  we  can  reduce  the  deter¬ 
mination  of  the  wave  functions  to  the  resolution  of  one  system  of  two  homogeneous 
coupled  integral  equations.  The  splitting  of  the  wave  matrix  into  S  and  D  parts| 
becomes  trivial. 

In  the  third  section,  we  use  various  methods  to  solve  these  equations. 

The  eigenvalue  of  the  integral  kernel  possesses  only  discrete  values  and  the 

lower  one  is  proportional  to  the  coupling  constant  of  nuclear  forces.  The 

computations  performed  in  the  ladder  approximation  leads  to  values  of  •£—  less 

4n 

than  8  (experimental  value  ~  15).  The  inclusion  in  the  interaction  kernel  of 
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fourth-order  terms  appears  as  necessary  a3  prohibitive  from  a  practical  point  of 
view. 


In  the  last  section,  we  present  a  simplified  method  to  obtain  a  relati¬ 
vistic  wave  function.  The  relativistic  corrections  are  divided  into  two  types. 

The  kinematical  corrections  and  the  properties  of  symmetry  correctly  describe 
the  spin  structure  of  the  Deuteron  which  is  represented  ty  an  entirely  known 
matrix.  The  relativistic  corrections  about  the  dynamics  of  the  system  are 
partially  included  In  a  scalar  function  which  is  equivalent,  in  a  non  relativistic 
limit,  to  the  Hul  then's  wave  function.  We  hope  this  solution,  essentially  pheno¬ 
menological,  is  reasonable  and  sufficiently  simple  to  permit  calculation.  The 
corresponding  percentage  of  D  state  in  the  deuteron  is  4  °/o  in  agreement  with 
experimental  results. 


II  INTEGRAL  ECUATIOJ3  for  the  DEUTERON  WAVE  FUNCTION 


II.  1. 


System  of  two  bound  nucleons  with  spin  ^ 


The  nucleon  is  represented  by  a  four-component  spinor  and  the  two- nucleon 
wave  function0(p)  i3  a  4  x  4  matrix,  the  'direct  product*  in  a  loose  sense  of 
the  two  spinors  fcr  the  two  particles. 


This  matrix  satisfies  the  Bethe  and  Salpcter  equation 


(j)  (p)  =  “ 


(y(l)Pl  +  iM)(Y(2)p.  +  iM)  m  f 

~n  -  Y  -  yV  / W(p,p» )  (|)  (p1)  dp« 

(Px  +  M  )(p2  +  M  ) 


(1) 


if  we  consider  a  PS  (PS)  coupling. 


The 


four-moments  p. ,  p„  are  related  to  the  moment  of  the  center  of 

AAr  1 


mass  P  and  the  relative  moment  p  through  the  equations  i 

aV\  *aa 


p  -  Pi  +  p2 

/w  vv  w* 


VVN  VAA  W* 


Po  = 


or 


2 


p 

--P 

2  w 
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The  matrices  y  ^  relative  to  the  particle  (i) 
i  =  1  and,  after  transposition,  to  the  right  if 


act  on  (b(p 

i  =  2  »>. 


(p)  from  the  left  if 


Y^P)  -  Yn 


Y 


y  being  an  usual  Dirac  matrix  and  y  f  its  transpose* 
(i  p* 


Let  us  introduce  the  B  matrix  : 


T  -1 
y =  B  V  B 


and  the  wave  function 


U’(p)  =  <b(p)  b"1  , 

*  1  vSA 

then  we  get  the  simple  relations  s 


Y(;V 

r'JV  tv ' 


The  new  wave  matrix  "^(p)  satisfies  the  integral  equation  : 

1  V\A 


y  (?)  = 

i*  /VSA 


a 


<»*  p2  *2 


W(p,p')  J"('yp1  +  iM)  y5'4i(p')  Y5WP2  +  ^  d£* 

(2) 


Let  us  write  YLKp) 

1  A r* 


in  a  reduced  form  having  four  elements  which  are  2X2  matrices 

/Yu 

'^21  ^22^ 


4)  Louis  de  Broglie  s  Thdorie  gdn4rale  des  parti cules  k  Spin  (M4thode  de  fusion) 
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Let  us  now  express  the  2x2  matrices  in  terms  of  the  unit  matrix  and  the  Pauli 
matrices 


^  (p)  =  S  (p)  I  +C?V  (p) 

-LJ.  y\  \  ±  VW  vs-~ 

VJ  (p)  =  S  (p)  I  +  vl(p) 

\fV\  *-  ws  C  X'A 


(p) 


u 

J21 


(p) 


b(p)  +  c(P) 


2 

B(p)  -  C(p) 
2 


-  -ft  )  ^ 

_*  F(p)  +  G(p) 

I  +o  — - — 

2 

_»"?(p)  -  G(p) 

I  +cr — ^ - — 

2 


If  we  substitute  these  expressions  for  the  reduced  elements  of  'Xl'(p)  in  equ,  (2). 


A** 


we  get  a  system  of  coupled  homogeneous  integral  equations  which  has  already  been 
written  down  for  the  scattering  problem 


Using  the  same  method  as  for  the  problem  of  two  particles  without  spin'  ' 
we  expand  the  scalar  functions  S  (p),  S  (p),  B  (p),  C  (p)  in  hyperspherical 

1  C.  >W  rvAA 

harmonics  and  the  vector  functions  in  vectorial  hyperspherical  harmonics.  We 
assume,  as  was  done  before,  that  the  interaction  kernel  W  (p,p!)  depends  only 
on  the  lengths  of  the  four  vectors  p,  p1  and  on  the  angle  (h)  between  them. 


This  assumption  allows  us  to  integrate  over  the  angles  ©  and  ,  the 
hyperspherical  harmonics,  both  scalar  and  vectorial,  being  orthonormalized. 

When  performing  the  integration  over  the  third  angle  p  we  are  led  to  introduce 
a  lot  of  auxiliary  functions  which  are  the  elements  of  a  16  X  16  matrix  which 
we  call  K(p  ;  J,  n,  n!)*  This  matrix  will  eventually  act  on  the  16  components 
vector  ^(p)  which  was  written  as  a  4  x  4  matrix  itself.  In  the  case  of  nucleons 
without  spin,  K  was  shown  to  be  simply  the  function  E^(p  \  1,  n,  n*)» 

We  order  the  16  components  of  ^(p)  in  an  arbitraiy  way  ;  the  non  vanishing 
elements  of  K  indicate  between  what  radial  components  of  \P(p)  a  coupling  does 

(3) 

exist.  We  have  established  in  a  former  work  on  the  nnol eon-nucleon  scattering 


the  following  results  : 
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a)  The  K  matrix  is  expressed  in  terms  of  four  8x8  sub-matrices.  The 
two  non-diagonal  ones  are  zero,  the  two  diagonal  ones  refer  to  a  definite  value 
of  the  parity,  the  total  angular  momentum  J  being  a  constant  of  the  motion. 

b)  The  sub-matrix  of  parity  (-1)^  i3  in  turn  expressed  in  the  same  w e~j, 
in  terns  of  two  4x4  sub-matrices  lying  along  the  diagonal  referring  respective¬ 
ly  to  singlet  (J  =  L,  m  =  0)  and  triplet  (j  =  L,  m  =  ±  l)  states,  m  is  the 
azimuthal  quantum  number. 

II,  2.  The  deuteron  case  j  J  =  1  state 


For  a  deuteron  in  the  ground  state,  J  =  1  and  parity  is  +  1.  The  wave 
function  will,  a  priori,  have  8  components  and  K  matrix  will  be  the  sub-matrix 
8  X  8  of  parity  +  1, 

Let  us  write  in  the  following  order  the  radial  components  A(p  ;  J  L  m  n) 
of  ^(p),  omitting  the  index  J  equal  to  1  and  the  index  L  when  L  is  equal  to  1 


Vx(p  ;  2  m  n) 
B  (p  ;  m  n) 


V^(p  ;  0  m  n) 
G  (p  ;  m  n) 


;  V2(p  ;  2  m  n)  ;  V2(p  ;  0  m  n) 
;  C  (p  }  m  n)  F  (p  ;  m  n). 


(?) 


Cn  the  other  hand,  the  matrix  K(p  ;  n  n')  is  an  integral  over  the  angle  (3  cf 
another  matrix  R(p,  (3)  : 

rn  o  n 


KX)i(p  ;n,  n')  = 


(p)  a^(p»  P)  ^'(P) 

0  (p2  +  M24o2-p2r2cos2|3 


sin^  p  dp 


(4) 


The  indices  i?'  are  related  respectively  to  the  components  X,  p..  With  the 
above  choice  (5)  for  numbering  the  radial  components  A^(p  ;  1,  m,  n)  we  have 

l  =  0  for  X  =  2,  4 

1  for  X  =  5,  6,  7,  8  (5) 

1=2  for  X  =  1,  3  . 


The  matrix  R  is  given  in  the  table  I  . 


Table  I  ;  R  Matrix  -  6  - 


sll-- 

2  2 
p  sin  p 

3 

R12  ~  2  ^2  Rn 

r 

™  ~  P  cos  p  +  (2 H  —  ]$] 

'1  R14“° 

v  =  \/ 

|  p  sin  p  (2M  -  |) 

R.,  =  -2-  R 

16  ^  15 

=  i  v|  P2  sin  P  COS  P 

r  =  4=  r37 

is  17 

R21  R12 

R22  =  "  R11 

R  =  0 

23 

R24  =  R13 

R25  =  Rlo 

R26  “  ~  R15 

R27  =  R18 

R28  =  "  B17 

Si  "  ' 

2  ?..  B2 

P  COS  P  4  - 

4  J 

^  =  0 

*33  "  R11 

^4  =  R12 

*35  H 

/|  P  cin  P  \ 

*37  =  R28 

R58  =  '  R27 

R41=0 

R42  -  Si 

R43  =  R21 

R44  =  +  R22 

R45 

*35 

R46  =  S$ 

R47  =  Se 

R48  =  R17 

*51  =  " 

2  *35 

*52  "  2  ^6 

S3-2  *15 

^4  "  "  2  R16 

%5  =  p2 

2 

—  V 

1 

^6=° 

=  -  2  i  p  M  cos  p 

^8  =  0 

& 

O ^ 

it 

4* 

2 

R62  =  “  R51 

R63  =  “  ^4 

R84  "  R55 

R65  =  0 

R66  =  “  ^5 

R67  =  ° 

RC3  =  ~  Si 

V- 

2R17 

P72  =  “  2  R18 

*13"  2  R17 

*74  =  2  R18 

4* 

\Jl 

1! 

I 

4* 

-a 

**76  ~  0 

R77  =  -  P2  cos2p  4  y2 

*78=° 

R81  “  "  2  R1S 

^2  =  2  R28 

Rq5  =  -  2  R1Q 

^  =  2  R17 

R86  =  *57 

*B8  =  ~  *71 

-  7  - 

The  selection  rules  due  to  the  parity  of  the  Gegenbouer  polynomials  show  that 
certain  elements  of  K  will  be  zero  after  integration  over  p.  Thus,  the  com¬ 
ponents  A^  Ag,  Aj,  A^,  A^t  Ag  will  be  real  while  the  components  A^,  Ag 
will  be  imaginary  j  we  have  n  +  n'  =  i  +  £ '  +  2  q  for  the  former  and 
n  +  n'  =»  ■£+  £'  +  2  q  +  1  for  the  latter. 

Then,  we  obtain  for  the  radial  components  A^  (p  ;  m,  n)  an  infinite 
system  of  integral  equations  coupled  by  the  indices  n  and  X  , 

_  _  roo 

An  (p  j  m  n)  =  v.  (p  ;  n  n')  /  An,(p,P*)  A  (p*  }  m  n')  p'5  dp' 

n*  n  H  J0  P 

(6) 


III  APHtOXIMATE  RESOLUTION  of  THB  SYSTEII  OF  INTEGRAL  EQUATIONS 

Using  the  same  approximation  methods  as  in  the  preceding  article  ^ 
we  propose  to  find  a  solution  to  the  system  of  integral  equations  (6).  We  saw 
that  these  equations  we’-e  very  weakly  coupled  with  respect  to  the  index  n  and 
that  the  wave  functions  fell  off  very  fast  as  n  increased.  The  function  W(j> ,gj) 
which  is  the  scalar  part  of  the  interaction  without  the  Dirac  matrices  corres¬ 
ponds  to  the  same  function  W(jj>,j3')  introduced  for  the  simplified  scalar  problem. 
Thus  it  turns  out  that  the  functions  An  ( P» p* )  are  the  same  in  the  two  cases. 
Actually  we  have  employed  the  ladder  approximation  but  the  arguments  for  the 
convergence  of  series  with  the  index  n  remain  valid  for  physically  reasonable 
interactions.  The  falling  off  of  An  (p»P')  with  n  is  very  fast,  we  shall 
decouple  the  equations  and  conserve  for  a  given  component  only  the  term  of  the 
smallest  order  in  n. 

III.l.  Terms  n  =  0  and  n'  =  0 

The  relation  connecting  ,  ft’,  n,  n'  and  q  shows  that,  in  this  case, 
only  two  components  intervene  :  V^(p  ;  0  0  0)  and  Vg(p  ;  0  0  0)  which  we  call 
V~  (p)  and  V~  (p).  The  system  of  integral  equations  (6)  is  than  reduced  to  the 
two  coupled  equations  t 


*■  8  — 


V“  (p)  =  a 


V:  (p)  =  a 


K22  (p) 


k42  (P) 


roo 


'  oo 


A0(p,p')  Vj(p')  p'3dp'  +K24(p)|  Aq(p,p')  V“(p!)  p'3dp’ 

0  •'0 

f<*>  /'OO 

A0(p*p')  v~(p')  p'pdp'  +K  (p)  A0(p,p')  v-(p’)  p'^dp* 


Jo 


(7) 


Henceforth,  by  convention,  K^^(p)  is  the  function  K  ;  n,  n1) 

corresponding  to  minimum  values  of  n  and  nf  for  which  it  is  not  zero.  Their 
calculation  becomes  that  of  the  functions  E^(p  ;  t  0  0). 

K22<p>  -  k44(p)  -  A1  -  y>  E0 


Vp)  *  - 


K42(p)  =  - 


-2-  (l  +  h(p))  +  4  M2 
4 

2  /  2  v  2 

-  (l  +  h(p))  +  UL  \ 

4  \4M  / 


E0(p) 


E0(P)  • 


Recall 


eJp)  = 


(p2  +  Y2)  [i?2  +  Y2)  +4p2M2^ 


1/2 


y  j  f  p  n  j 

2  o  .  1  l/2 


h  (?)  = 


(p2  +  72)  -  {(p2  +  y2)  +  4  p2  M2  \ 
(p2  +  Y2)  +  {(p2  +  y  )  +  4  P2  I M2} 


1/2 


and  the  binding  energy  of  deuteron  -B  ippears  tlirough  the  parameter  y  thus 
defined 


2 

Y 


=  MB. 


To  facilitate  calculation  we  define  two  auxiliary  functions  I^(p)  and  l2(p) 


*1,2  (p> 


roO 

A0(p,p’)  v“>2(P>)  p'3 


Jo 


dp' 


(3) 


permitting  us  to  write  the  system  (7)  as 


(9) 
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v‘  (p)  =  <x^k22(p)  ix(p)  +  k24(p)  i2(p)J 
Vj  (p)  =  ajK42(p)  I^P)  +  K22(p)  I2(p)| 

III.  2,  Terms  n  =  1,  n'  =  0 


At  this  order  of  the  approximation,  two  now  components  of  the  wave 
matrix  B(p  ;  1  u  l)  and  G(p  ;  1  m  l)  appear.  We  call  them  B  (p)  and  G  (p)  and 
express  them  with  the  help  of  1^  and  I2  . 


B  (p)  =  a  |k52(p)  Ix(p)  +  K^p)  I2(p)j 
G  (p)  =  a  |k62(P)  Ix(p)  +  K64(p)  I2(p)| 


(10) 


The  four  new  functions  K^^(p)  are  then  given  in  terms  of  E^(p)  by 

n  \ 2 


vp>  =  y 2  p  ” Ei<p) 


62 
where 


h(P)  , 

\(v)  =  (1  -  —  )  Eq(p) 


«54(p)  =  -  2  P  M  El 
K64(p)  «  -  2V2  p  H  Ex(p) 


We  can  invert  the  equations  (9) 

11 = y?22  vi  ■  ** 

12  ■  y  [-  k  « ri  *  ^  y 

where  we  put 

W  (p)  =  K22  K44  -  K24  K42 
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and  see  that  B  (p)  and  G  (p)  are,  at  this  point  in  the  calculation,  simple  linear 
combination  of  V~  (p)  and  (p),  the  coefficients  being  expressed  in  terms  of 
the  functions  K^(p). 

Finally  one  should  notice  that  the  binding  energy  of  the  Deutercn  being 

/  Y  \2  1 

very  small  (3  =  2,23  MeV),  the  number  —  )  is  much  less  than  unity  (  ~  5 

\2M/  2  °'° 

K^-2(p)  and.  (p)  are  thus  negligible  before  K^(p)  and  K^(p).  Doing  this, 

we  find  a  simple  relation  between  B  and  G 

G  (p)  =  2 Vi  (p)  (ll) 

which  is  certainly  valid  to  a  very  good  approximation, 

III.  3.  Terns  n  =  0,  nf  =1 


The  equations  (7)  are  now  modified  by  the  addition  of  a  coupling  term 
between  V~  (p),  V~  (p)  and  B  (p),  G  (p).  We  thus  calculate  a  first  correction 
to  the  solutions  of  the  equations  (7). 


If  we  introduce  the  auxiliary  functions  I^(p)  and  I^(p) 


I-j(p)  =  (  A1(p,p')  B  (p')  p'3  dp' 


I4(p)  =  f  A^P,?')  G  (p')  p'3  dp'  , 

^0 

the  system  (9)  becomes 


V-  (p)  =  * 
Vi  (p)  =  a 


K22  I1  +  K24  I2  +  K25  H  +  K26  I4 


K.„  I.  +  K  I„  +  K.c  I,  +  K,.  I, 
42  1  44  2  45  3  46  4 


(12) 


with 


K25(p)  =  P  M  E1 


K45(p) 


=  ./lf  L(,) 


2M 


I  25 


K26(p)  =  -\/2  K25(p) 
K46(p)  *  (U  K26<>> 
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The  problem  Is  then  to  solve  the  coupled  system  (10)  and  (12).  By  eliminating 
B  (p)  and  G  (p)  we  obtain  two  coupled  integral  equations  for  V”  (p)  and  V2  (p) 
alone.  The  difference  between  this  system  and  (8)  will  enable  us  to  test  our 
development  in  series  with  respect  to  n. 


But,  it  it  legitimate  to  neglectQ 
then  valid  and, the  system  (12)  is  equivalent  to 


before  unity. 


Relation  (ll)  is 


ri  +  *24  - 


25  "5 


(13) 


The  seoond  equation  (9)  is  not  modified,  so  that,  taking  (8),(9),(ll) 
into  account,  we  get  only  one  substitution. 


k24(p)  A0(p,p')  c=>  k24(p)  A0(p.p') 

roc  - 

-  ^(p)  A^p.p”)  K^p")  ^0(P",  P')  P"^  dp’» 

J0 

which  can  be  easily  calculated  in  the  plane  (p,  p')  since  all  the  elements  are 
known. 


III.  4»  Terms  n  =  2,  n'  =0  (real) 


In  this  case ,  the  orbital  angular  momentum  components  of  order  2  appear. 

They  are  V, (p  ;  2  m  2)  and  V.(p  ;  2  m  2),  which  we  call  V,+(p)  and  V«+(p) 
12  12 

respectively.  They  are  easily  known  once  V7  (p)  and  V~  (p)  have  been  calculated. 


V1+(p)  =  a  K12(p)  I1(p) 
V2+(p)  -  a  K^p)  I2(p)  . 


(14) 


The  two  functions  K^(p)  and  K^4(p)  are  equal 

\f 5  ? 

Ki2(p)  =  K-j^p)  mm~T  p 


h  (p) 
2 


h2(p) 


10 


Eo  (p) 
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IV 


III,  5.  Terms  n  =  2  n'  =0  (imaginary) 


The  two  imaginary  functions  that  we  introduce  here  C(p  ;  1  m  2)  and 

2 

F(p  ;  1  m  2)  are  to  be  multiplied  by  the  function  ^  ^  (p).  It  is  easy  to  see 
that  in  configuration  space,  the  same  thing  happens.  At  the  limit  of  equal  times, 
P  =  2  and  (*[)  =  0.  For  that  reason,  we  do  not  worry  about  these  imaginary 
functions. 

III.  6.  Conclusion 


Taking  into  account  in  each  case  only  the  lowest  terms  in  n  and  n* ,  we 
are  led  to  the  following  programme 

a)  Solve  the  system  (7)  to  obtain  V~  (p),  V~  (p),  1^  and  , 

b)  Calculate  B  (p)  and  G  (p)  by  equations  (10)  , 

c)  Calculate  V^+(p)  and  V2+(p)  by  equations  (14)  . 


SEPARATION  into  STATES  and 


The  wave  matrix  ill  (p)  in  our  particular  case  has  therefore  the  form  : 
G  .  V  (p) 

i  'W 


'll)  (?)  = 


B  (jP,)  aG  (p) 

- L.  +  - — 


2 


B  (p)  O'G  (p) 


G  .  V  (p) 
c  **• 


In  the  state  of  total  angular  momentum  J  =  1,  it  becomes  : 


m  0 


m  0 


Vl,2  “  71,2  ^  "101  ’  M  +  vi,2  ^  “(  2  1  (e'f'  W 

b  (p)  -  b  (p)  (p)  y^e.f) 

— *  « 

-»  c\J  rn  0 

g  (p)  -  o  (p)  .  ^lxl  (e,<p,  P) 
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V  m  n 

The  vectorial  hyper spherical  harmonics  J  (0fip,  P)  ^y  be  obtained  from 

J  Jj  JL  ' 

scalar  hyperspherical  harmonics  by  the  relations  : 


y  m  n  _  i  ?  v  m  n 

JJ  J  r  \rz7T-rri  L  ,yJ 


Vj(j  +  1) 


y 


mfn 


j,  J+i,  i  y'(j  +  i) ( 2  j  +*i) 


— » 


(J  +  l)  e  +  i  e  X  L 
P  P 


% 


m,n 


V 


m5n 


J,  .>-1,  1  y|j  +  j  + 


Je^-ie^XL 


ym>n 


We  easily  deduce  the  two  relations 


y 


m  n 


(o!  T)  .  7j 


i/jLu  ym'G  +  \/-  j  ym,n 

V  2  J  +  1  "yJ,  J+l,  1  V  2  J  +  1  'yJ,  J-l, 


y,  ‘  1  -  (?  !p  •  *  l/rrfi  V,"  “+1, !  -  Vj”  lh  ! 


using  the  property  e  L  =  0,  e  being  the  unitary  vector  in  the  direction 'pi 
P  P 


Applying  these  equations  for  J  =  lf  the  wave  matrix  becomes  simply  : 


^  (?)  =  Mg(p.  P)  <?  yx  q  ^e.vp)  +  m^p,  P)  (?  Vx  2  ^e,^)  . 


Here 


V. 


m 


J  L  1 


are  ordinary  vector  orthonormal  spherical  harmonics. 


The  wave 


matrices  for  the  S  and  D  states  are  given  by 
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=  \fl 


v"  (p) 


j[b(P)-\|2  G(p)jsinpa< 


\  [b(P)  +  V2  G(p)J  sin  p  ae 

— pr  sin2p  V  +(p) 

VI  1 

\  |j(p)  -  G(p)j  sin  Pa 


V'(P) 


j  B(p)  +  G(p) j  sin  p  a 


~  sin2p  V2+(p) 

\f5  * 


It  is  then  easy  to  calculate  the  proportion  of  S  and  D  states  by  the  integrals 


Trace J  |Mg|2  P^  dp  sin2p  dp 
Trace  f| ^(p) | 2  dp 

J  ***  VNA 


Trace  |  I I ^  dp  sin2p  dp 


PD  = 


Trace 


The  integration  over  the  angle  p  are  immediate.  After  calculating 


the  traces  we  obtain 


ns  +  nd 


ns  +  nd 


n.- 


|v-(p)|2+  |v"  (p)|2  +  l  |b(p)|2+  ^  |g(p^  p3  dp 


Vl+(p)|2  +  |V2+(p)|2  +  ?  |B(p)|2+  5  |G(pf  P5  dP 
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V  NUMERICAL  RESOLUTION  of  INTEGRAL  EQUATIONS  in  the  LADDER  APPROXIMATION 

As  we  saw  before,  the  possibility  to  decouple  the  equations  with  respect 
to  the  index  n  focuses  all  the  difficulties  on  the  solving  of  the  system  of 
coupled  integral  equation  (7).  After  fixing  the  binding  energy  of  deuteron  at 
its  experimental  value,  we  intend  to  determine  the  coupling  constant  and  the 
wave  functions  (p)  and  V~  (p)  from  which  we  can  express  all  the  others  by 
quadratures. 


V.  1.  Approximate  resolution  : 


We  substitute  for  the  exact  kernel  An(p>P')  an  approximate  and 


separable  one  : 


zyP,p')  <=>  — 


choosing  for  q  the  average  momentum  of  nucleons  in  the  deuteron  (qc^i  p).  2 

^  £r 

We  can  then  solve  the  new  system  exactly  and  we  find  a  =  0,13  ieS,  _f 

the  coupling  constant  of  nuclear  forces,  is  equal  to  16. 


V.  2„  Kellogg  method  of  iteration 


We  have  made  a  numerical  calculation  on  IBM  650  to  solve  the  equations 
by  successive  iterations.  We  thus  define  a  series  of  values  for  a  at  each  stage 
of  the  iteration.  The  result  is 


a  ^0,046 


i.e. 


-  —  6  , 


and  the  proportion  of  D  state  is 


1  /o . 


We  have  verified  that  the  development  with  respect  to  n  is  excellent 
as  well.  The  solving  of  the  system  (lj)  instead  of  (?)  leads  us  to 


<XCZ  0,048 
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i.e,  ^  ^ 6 ,5  • 

4n 

We  then  obtain  a  clear  disagreement  with  the  experimental  value 
deduced  for  example  from  meson-nucleon  scattering. 


2 

4 ft 
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The  most  reasonable  explanation  for  this  seems  to  be  the  inability  of 
the  ladder  approximation  to  describe  this  phenomenon.  It  would  be  necessary  in 
a  meson  theory  to  include  fourth  order  teims,  Unfortunately  this  calculation 
is  inextricable  and  we  shall  not  do  it  although  the  formalism  remains  valid. 


VI  APPROXIMATE  WAVE  FUNCTIONS 


The  calculations  with  a  kernel  including  the  fourth  order  term  being 

inextricable,  we  have  modified  our  method  in  order  to  obtain  a  deuteron  wave 

(5) 

matrix  available  in  concrete  problems  , 


If  we  consider  Bethe  and  Salpeter's  equation  : 


a 


(•/ -  iM)(-/2)p2  -  iM) 


•'Z1)  Y^  f w  (p,p’)  CD  (p* )  dp'  , 

5  5  J  V.  w>  vW  VvA 


we  can  see  that  there  are  ?  types  of  relativistic  corrections  : 


-  the  kLnematical  corrections  take  the  relativistic  energy  and  the  spins 
of  the  nucleons  into  account,  they  are  included  in  the  expressions  (y^p^  -  iMj 
and  (y^p2  -  iM)  ; 


,(i) 


-  the  properties  cf  symmetry  are  correctly  represented  by  the  matrices 


"Y  _  which  show  that  we  have  used  a  pseudo-scalar  mesor.-nucleon  interaction  ; 
5 


-  the  dynamical  corrections  are  included  in  the  expression  of  the  kernel 


W  (p,p'). 


In  our  previous  calculation  ,  we  have  precisely  taken  the  first  two 


5)  We  are  indebted  to  Professor  M.  L£vy  for  this  suggestion. 
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kinds  of  correction  into  account.  On  the  other  hand,  we  have  only  partially 
considered  the  relativistic  effects  in  the  dynamical  corrections,  assuming  the 
hypothesis  that  the  interaction  depends  only  upon  the  square  of  the  space-time 
distance  between  the  two  particles.  In  particular,  we  have  used  the  ladder 
approximation. 

In  separating  the  two  types  of  corrections,  we  have  obtained  : 

a)  a  correct  structure  of  the  wave  matrix  -  i.e.  the  distribution 
between  it 3  different  components  -  due  only  to  the  spin  of  the  nucleons  and  to 
the  pseudo-scalar  character  of  the  interaction, 

b)  the  wave  functions  (p),  (p),  which  are  approximative,  but 

which  take  some  relativistic  effects  into  account. 


In  order  to  simplify  the  integral  equations,  we  keep  only  the  leading 
terns ,  and  this  enables  us  to  obtain  a  simpler  structure  of  the  wave  matrix. 
Considering  that 


V  \2 

3)  «l 


and 


V"  (p)  V-  (p)  =  -  b2(p) 


with  a  very  good  approximation,  we  can  obtain  the  wave  matrices  and 


“s 


4 


M  =-4- 

D  v* 


§5  (l  -  5)  sin  p  a  T 


-  -  V 


V~  (p) 


1  -  -  2  2 
8  2  /.  hx  , P  v  ,  2. 

3  U  "  3)  '2M'  S^n  ^ 

?  1  +  h  p 


|  (|)  (l  -  |)  sin  P  O'  ej 


|  u-!>  (9)  ^2p 


Vj(p) 


(18) 
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The  only  remaining  unknown  scalar  function  V”  (p)  corresponds,  in  the 
non  relativistic  approximation,  to  the  wave  function  for  a  deuteron  without 
spin.  We  represent  it  graphically  in  fig.  I. 

The  calculation  of  the  immixture  of  D  state  gives  4  °/o,  which  is  in 
satisfactory  agreement  with  experiments.  We  think  that  this  very  phenomeno¬ 
logical  determination  of  the  deuteron  wave  function  is  better  the  preced¬ 

ing  one.  It  amounts  to  writing  the  wave  matrix  with  the  help  of  a  scalar 
function  resembling  that  of  Hulthen,  and  known  matrices  describing  a  structure 

bound  to  the  existence  of  nucleon  spin.  In  other  words,  the  matrices  M  and  M. 

3  3  S 

express  the  coupling  between  the  states  and  D^  of  Deuteron  and  include 
relativistic  effects  inside  the  Deuteron  which  are  included  in  (p),  a  func¬ 
tion  of  the  same  type  as  (p),  the  solution  corresponding  to  the  problem 
without  spins. 


VII  CONCLUSION 


We  have  developed  a  formalism  to  determine  the  relativistic  deuteron 
wave  function  and  the  coupling  constant  of  nuclear  forces  from  the  binding 
energy  of  the  deuteron  which  has  been  experimentally  determined. 

We  have  used  the  interaction  given  by  the  meson  theory  about  nuclear 
forces  in  the  ladder  approximation  of  second  order.  The  theoretical  results 
make  us  foresee  a  coupling  constant  twice  too  small  and  this  does  not  allow  us 
to  take  too  seriously  the  results  obtained  about  the  wave  function. 

We  have  tried  to  consider  the  mass  of  meson  n  as  a  phenomenological 
parameter.  The  most  favourable  results,  corresponding  to  p  =  0,1  (g^/4it~7 
and  p^Ci  1,5  °/o)  show  that  the  coupling  constant  as  well  as  the  proportion  of 
D  state  are  parameters  which  vary  extremely  slowly  with  p,  Besides,  this  has 
been  verified  in  the  diffusion  problem. 

We  have  also  thought  of  adjoining  an  arbitrary  phenomenological,  term 
respecting  the  conditions  of  invariance  to  the  mesonic  tern  of  second  order  j 
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an  attempt  has  been  made  but  was  not  very  conclusive  5  perhaps,  there  are,  after 
all,  some  attempts  to  be  made  in  this  direction,  but  the  justification  and  the 
theoretical  origin  of  such  terms  seem  to  be  extremely  problematical  and  finally 
doomed  to  failure. 

With  extremely  weak  approximations,  it  is  possible  to  give  to  the  wave 
matrix  a  structure  entirely  known  from  one  scalar  function  only.  This  method 
leads  us  to  a  very  simple  form  and  will  permit  numerical  computations.  Parti¬ 
cularly,  the  most  immediate  and  interesting  application  is  the  calculation  of 
the  elastic  electron  deuteron  scattering  cross  section  in  the  impulse  approxi¬ 
mation.  The  problem  becomes  very  clear  to  understand.  We  know  how  the  inter¬ 
action  acts  on  the  neutron  and  on  the  proton  and  all  the  trace  calculations  can 
be  performed  before  integration.  The  modifications  due  to  the  spin  of  the 
nucleons  will  appear  naturally  and  the  relativistic  corrections  will  be  rapidly 
evaluated. 


